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Abstract 

We derive Painleve-type expressions for the distribution of the m"* 
largest eigenvalue in the Gaussian Orthogonal and Symplectic Ensembles 
in the edge scaling limit. The work of Johnstone and Soshnikov (see |7|, 
|1UI '1 implies the immediate relevance of our formulas for the m*'' largest 
eigenvalue of the appropriate Wishart distribution. 



1 Introduction 

The Gaussian /3-enscmblcs are probability spaces on n-tuples of random vari- 
ables {Ai, . . . , Ajv}, with joint density functions 



Pn f}{^i, ■ ■ ■ , Ajv) — Pn — Cn/s cxp 



1 ^ 



niA,-Afci^. (1.1) 



The Cnp are normalization constants, and by setting [3 = 1,2,4 we recover the 
Gaussian Orthogonal Ensemble (GOEat), Gaussian Unitary Ensemble (GUEjv), 
and Gaussian Symplectic Ensemble (GSE^r), respectively. We restrict ourselves 
to those three cases in this paper, and refer the reader to |3] for recent results 
on the general (3 case. Originally the Xj are eigenvalues of randomly chosen 
matrices from corresponding matrix ensembles, so we will henceforth refer to 
them as eigenvalues. With the eigenvalues ordered so that > Aj+i, define 



/2N 



2-1/2 7V-1/6' 



(1.2) 



to be the rescaled m eigenvalue measured from edge of spectrum. A standard 
result of Random Matrix Theory about the distribution of the largest eigenvalue 
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in the /3-ensembles is that 

whose law is given by the Tracy- Widom distributions. 
Theorem 1.1 (Tracy, Widom [131. [T4|). 



(1.3) 



F2{s) : 



Pgc/b(-^i <s) = exp 



(a; — s) q^{x)dx 



.(Ai < s) 



Fii^) [P,,,(Ai < .) 



F2 ■ exp 



= F2 ■ cosli 



q(x)dx 



q{x)dx 



(1.4) 

(1.5) 
(1.6) 



The function q is the unique (see |nij[21) sohition to the Painlcve II equation 



q" =xq + 2q\ 



(1.7) 



such that q{x) ^ Ai(2;) as x 00, where Ai(2;) is the solution to the Airy 
equation which decays like ^ tt^^/^ x^^/* exp (— | x'^/^) at +00. The density 
functions fp corresponding to the are graphed in Figure 




Figure 1: Tracy- Widom Density Functions 

Let F2{s,m) denote the distribution for the to*'* largest eigenvalue in GUE. 
Tracy and Widom showed in J3] that if we define -^2(5, 0) = 0, then 



F2(s,m+1)-F2(s,m) 



(-ir 



■D2{s,X)\^^^, m>0, (1. 



*The square root of 2 in the argument of F4 reflects a normalization chosen in 11.11 to 
agree with Mehta's original one. It can be removed by choosing a different normalization. 
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where 



1)2(5, A) = exp 



(a; — s) q^{x, X)dx 



(1.9) 



and q{x, A) is the solution to H1.7|l such that q{x, A) ^ VX Ai{x) as x 00. An 
intermediate step leading to (|1.9() is to first show that D2{s, A) can be expressed 
as a Fredholm determinant 

i?2(s,A) =det(/- A/CAi), (1.10) 
where /Cai is the integral operator with kernel 

x-y 

In the P = 1,4 cases a result similar to (|l.lUf) holds with the difference that the 
operators in Dp{s,X) have matrix-valued kernels (see e.g. JJ). In fact, the 
same combinatorial argument used to obtain the recurrence 1)1. 8|l in the (3 = 2 
case also works for the /3 = 1, 4 cases, leading to 

Ff^is, m + 1) - Ff^is, m) = d'/\s, A)|^^^ , m > 0, /3 = 1, 4, 

(1.12) 

where i^^(s,0) = 0. Given the similarity in the arguments up to this point 
and comparing (|1.9() to (|1.4|l . it is natural to conjecture that _D^(s,A),/3 = 
1,4, can be obtained simply by replacing q{x) by q{x,X) in (|1.5(l and l|1.6(l . 
However the following conjecture, which had long been in the literature, and 
whose verification is the content of Corollary 1)2. 2|l . hints that this cannot be 
the case: 

Conjecture 1.2 (Baik, Rains ^). In the appropriate scaling limit, the distri- 
bution of the largest eigenvalue in GSE corresponds to that of the second largest 
in GOE. More generally, the joint distribution of every second eigenvalue in the 
GOE coincides with the joint distribution of all the eigenvalues in the GSE, with 
an appropriate number of eigenvalues. 

Forrester and Rains subsequently proved (see jS]) the equivalence of alternate 
GOE eigenvalues and GSE eigenvalues at finite N ensemble level, lending weight 
to Conjecture \1.2\ . This so-called "interlacing property" between GOE and 
GSE had been noticed by Mehta and Dyson (see |H]). Conjecture H1.2|l does not 
agree with the formulae we postulated for Dp{s, X), (3 = 1,4. Indeed, combining 
the two leads to incorrect relationships between derivatives of q{x, A) evaluated 
at A = 1. To be precise, the conjecture is true for 1)4(3, A) but it is false 
for Di{s,X). The correct forms for both Dp{s,X),l3 =1,4 are given below in 
Theorem (|TT|) . 

This work also extends that of Johnstone in [7] (see also @]), since Fi{s,m) 
gives the asymptotic behavior of the m*'* largest eigenvalue of a p variate 
Wishart distribution on n degrees of freedom with identity covariance. This 
holds under very general conditions on the underlying distribution of matrix 
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entries by Soshnikov's universality theorem (see ^Hl for ^ precise statement). 
In Table ^ we compare our distributions to finite n and p empirical Wishart 
distributions as in [7|. 



2 Statement of the Main Results 

Theorem 2.1. The distributions for the m*'* largest eigenvalues in the GOE 
and GSE satisfy the recurrence with 

A-l-cosh^(s,A) + yisinhpt(s,A) 
Di{s,\) ^ D2{s,\) ^— ^ , (2.1) 

Di{s, A) - D2{s, A) cosh^ f , (2.2) 



where 

/oo 
q{x,X)dx, A:=2A-A^ (2.3) 

and q{x, A) is the solution to (|1.7|l such that q{x, A) ~ \/A Ai(.T) as x — *■ oo. 
Corollary 2.2 (Interlacing property). 

F4(s,m) = Fi(s,2m), m > 1. (2.4) 



In the next section we outline the proof of these theorems. In the last, we 
present an efficient numerical scheme to compute F/}{s,m). We implemented 
this scheme using MATLAB, and compared the results to simulated Wishart 
distributions. 



3 Sketch of the Proofs 

3.1 Distribution for the Next Largest Eigenvalues and 
Finite Gaussian Ensembles 

With the joint density function defined as in Hl.l|) . let J be an interval on the 
real hue, and X = Xj (x) its characteristic function. We denote by X = 1 — X 
the characteristic function of the complement of J, and define X^ ~ 1 — AX. 
Furthermore, let Ep^Nijn, J) equal the probability that exactly the m largest 
eigenvalues of a matrix chosen at random from a (finite N) /?-ensemble lie in 
J. We also define 



X^{xi) ■ ■ -X^ixN) Pn f3{xi 



,xn) dxi ■ ■ ■ dxN- (3.1) 
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Figure 2: 10^ realizations of 10^ x 10'^ GOE matrices; the solid curves are, from 
right to left, the theoretical limiting densities for the first through fourth largest 
eigenvalue. 



For A = 1 this is just £'/3_n(0, J), the probability that no eigenvalues lie in J. 
The following propositions are easy combinatorial facts that can be proved by 
induction (see e.g. |12|). 

Proposition 3.1. 

^P^UJ) = J2^-^)'' (Vj f [ Pnp{xu ■■■,XN)dxi-- -dxN. (3.2) 
fe=0 \ ' / J J 

Proposition 3.2. 



(-1) 



771 Jrn 



TO > 0. 



(3.3) 



A=l 



The next step is to find a useful expression for the multiple integral Eg^l^(J). 



It turns out that through standard RMT techniques (see e.g. Jl]), the inte- 
gral can be expressed as the determinant of an operator on the Hilbert space 
L2(J) X L^{J). Let 



for 



S + 'ip(g>e(p SD -ip^ip 



(3.4) 
(3.5) 
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Here e is the integral operator with kernel e{x — y) = | sgn(x — y), and D denotes 
the differentiation operator S is the integral operator with kernel 

x-y 



and the functions and -0 are 

'^(^) =\-w] Vn{x), (3.7) 



= f -j ^jv-i(x), (3.8) 

where 

^.{x)^-^^==e-^''-H,{x), (3.9) 

and the i?jv(a;) are the classical Hermite polynomials. This implies that the 
(Pn{x) are orthonormal with respect to the Lebesgue measure on M. Similarly, 
let 

Following the same approach as in jlSj and |14| . we arrive at 

= (3.11) 

3.2 Edge-Scaling 

3.2.1 The GOE case: reduction of the determinant 

The above determinants are Fredholm determinants of operators on L'^{J) x L^{J). 
Our first task will be to rewrite these determinants as those of operators on 
LP'iJ)- This part follows exactly the proof in ^3]. To begin, note that 

[S ,D]=Lp®i) + '4)®ip (3.12) 

so that (using the fact that D e ^ e D = I ) 

[£,5*] = £5*— 5£ 

= £5'L>£- £i:)S'£ = £ [S*,!?] £ 

= £ (g) £ + £?/' (g) (^£ 

= eif ® e*ip + eip ® (fi 

= —eip^eip — eip^eip, (3.13) 
where the last equality follows from the fact that £* ~ ~e. We thus have 
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D {eS + etp(E)eLp) = S + ip(E)e(p, 
D {eSD-e%l)®Lp) = SD-iP®(p. 

The expressions on the right side are the top entries of /Ci^tv- Thus the first row 
of /Ci.AT is, as a vector, 

D {eS + eip e(p,eS D - eip if) . 
Now imphes that 

e S + ei/j ^ eip = S e — eip ® eijj. 

Similarly H3.12|l gives 
so that 

tSD-eilj®ip^eDS + tip®i^ = S + tip®'4}. 
Using these expressions we can rewrite the first row of /Ci,7v as 

D (S e - e(p(g> eilj,S + e(p(g>xp) . 
Applying e to this expression shows the second row of JCi^n is given by 

Now use H3.13|l to show the second row of ICi^n is 

{S e - e + e(p(g)eil!,S + £(f(E)tp) . 

Therefore, 



XD 0\/ {Se-t(p®e^)X {S + etp®ip)X 
X ) \ {Se-e + eip®e'il))X {S + eip<g)^)X 

Since ICi^n is of the form A B, we can use the fact that det(/— ^ B) = det(/— B A) 
and deduce that Di^n{s, A) is unchanged if instead we take /Ci^at to be 



^1. 



N 



{Se-e(p®eip) X {S + e ip <S) ip) X \ f XD 
(Se-e + eipi^e^p) X {S + e (p ® ^j) X ) \ X 

[Se- e(p®ei;) XD (5 + e <^ (g) V) X " 
{Se-€ + eLp®ei))XD {S + e (p ® ifj) X 
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Therefore 



Di,n{s, A) = dot 



I - {Se-eip(g)eip) XXD - (S + e ip (g> Tp) XX \ 
-{Se-e + eifSieip) XXD I - {S + e (p <g} ijj) XX )' 

(3.14) 

Now we perform row and column operations on the matrix to simplify it, which 
do not change the Fredholm determinant. Justification of these operations is 
given in We start by subtracting row 1 from row 2 to get 

I - {S e- eip<^ e^p) XXD - {S + e ip <^ ip) XX 
~I + eXXD I 

Next, adding column 2 to column 1 yields 



I - {St~ ep^etp) XXD-{S + eip®ilj) XX -{S + etp<^i)) XX 
eXXD I 

Then right-multiply column 2 by —eXXD and add it to column 1 to get 



I - (Se-eifSeip) XX D + (S + e ip ® ijj) XX {eXXD - I) -{S + ep®ip) X 

/ 

Finally we multiply row 2hyS + eipi^ip and add it to row 1 to arrive at 



/ I-iSe-eif^eip) XXD + {S + eip®i;) XXieXXD-I) 
det 1^ / 

Thus the determinant we want equals the determinant of 



I - {Se- ep(g)eijj) XXD + {S + e(p(g)ilj) XX (eXXD - I) . (3.15) 

So we have reduced the problem from the computation of the Fredholm deter- 
minant of an operator on L^{J) x L'^{J), to that of an operator on L^{J)- 



3.2.2 The GOE Case: difTerential equations 

Next we want to write the operator in H3.15|l in the form 

(/ - IC2.N) (^I~Y.c^^® , (3.16) 

where the and (3i are functions in L'^(J). In other words, we want to rewrite 
the determinant for the GOE case as a finite dimensional perturbation of the 
corresponding GUE determinant. The Fredholm determinant of the product 
is then the product of the determinants. The limiting form for the GUE part 
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is already known, and we can just focus on finding a limiting form for the 
determinant of the finite dimensional piece. It is here that the proof must be 
modified from that in ^1]. A little simplification of (|3.15|l yields 

I - XSX- XS {1- XX) eXD - X {eip X^/j)- X {e^p XX) eXD. 

Writing e [ X , Z? ] + X for e X and simphfying (1 - A X) X to (1 - A) X gives 

I - XSX~ X {I - X) SX- X {eip ® Xijj) ~ X {1- X) {eip ® X^p) 

-XS (I-XX) e[X,D]-X{eip(g> ip) (I-XX) e[X,D] 
= I - {2X - X^) S X - {2X - X"^) {e(p (g) Xi^) - X S {I - XX) e [X , D] 
- A(ev? (g) ij) (1 - Ax)e [X,D]. 

Define A = 2 A — A^ and let ip ^ f, and ^ -0 so that XS ^ S and 

H3.15|l goes to 

I-SX-ieip (E)Xilj)~ts{l-XX)e[X,D] 
A 

-jieip® AX)e [X,D]. 

X 

Now we define i? ;= (/ — S X)^^ SX = (/ — SX)^^ — I (the resolvent operator 
of SX), whose kernel we denote by R{x,y), and Qe '■= {I — SX)~^eip. Then 
(|3.15|l factors into 

A={I-SX)B. 

where B is 

I-{Q, ® Xi:) - ^ (I + R) S {1 - XX) e [X , D] 
A 

- jiQe ® ^)il- XX)e[X,D], X^l. 
A 

Hence 

Di,jv(s,A) =i?2,Jv(s,A) det(S). 
In order to find det(i?) we use the identity 

2 m 

e [X,i5] =^(-l)^-efc(g)(Sfe, (3.17) 
fc=i 

where Cfc and S^. are the functions e(a; — aj.) and 5(.t — aj.) respectively, and the 
are the endpoints of the (disjoint) intervals considered. J ~ U^!L2(a2fc_i, a2fc)- 
We also make use of the fact that 

a(g) 6 • c(8) d = (6,c) • a® d (3.18) 
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where ( . , . ) is the usual L^-mncr product. Therefore 



2m 



(1- Ax)e [X,i?] =^(-l)'=Q.®iA-(l-Ax)efc(^4 

k=l 
2m 

= ^(-l)'=(V',(l-AX)efe) Q,® 4. 



fc=i 



It follows that 



equals the determinant of 



Di,n{sA) 
D2.n{s, A) 



I -Qe(E>Xip 



+ (1- AX)efe + (i^,(l- AX)efc) Qel^^. 

fc=i 

We now specialize to the case of one interval J ~ (t, oo). so m = 1. ai = i and 
02 = oo. We write et = ei, and eoo = £2, and similarly for 6k- Writing the terms 
in the summation and using the facts that 

eoo = (3.19) 

and 

(1 - AX) et = (1 - AX) + (1 - AX) X, (3.20) 

then yields 

/-Q, ^XV--^ [(^ + i?5)(l-AX) + (^,(l-AX)) Qe](^{St-5oo) 
2A 

+ ^ [(5 + (1 - AX)X + (i^ , (1 - AX)X) Q,] St 
A 

which, to simplify notation, we write as 

/ - g, ® XV- - ^ [(5 + i?^) (1 - AX) + ai,A Qe] «) {St - <5oo) 
2A 

+^ [(5* + i? 5) (1 - A X) X + ai^A Qc] ® 5* , 
A 

where 

ai,A = (^',(1- AX)), ai,A = (^'.(l- AX)X). (3.21) 
Now we can use the formula: 

det (^I-J2a,(g, p}j = det {Sjk ~ {a, , Pk ))i<j^k<L (3-22) 
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In this case, i = 3, and 



ai =Qe, a2^j [{S + RS){l-\X)+ai^xQe]^ 
A 

a^^~h{S + RS){l-\X)X + ai,A Q.] , 
A 

f3i = Xi>, /32 = <5t - <^oo, Ps^St. 
In order to simplify the notation, define 



(3.23) 



Q{x,X,t) := (/- 




P{x,X,t) := (/- 


-5X)-iV, 


Q,(x,A,t) := (/- 




P,(a-,A,<) := (/- 


-Sxy^etp, 


g„ :=g(t,A,t), 




p„ :=P(t,A,i), 




9, :==g,(t,A,t), 




~P,(i,A,i), 




:= (Q ,Xe(;9) = 




V, := (Q,XeV) 




i), := (P,Ae<^) = 






= (-Pe,XV), 



(3.25) 



Ql,A 



(1 - AX)Pda;, 
(1 - \X)Qdx, 



Ql,A 



(1 - AX)XPda;, 
(1 - AX)XQdx, 



Tli.x:^ j {I- \X)R{x,t)dx, ni.x.= j {I- \X)XR{x,t)dx. (3.26) 

Note that all quantities in H3.25|l and (|3.26|) are functions of t alone. Further- 
more, let 



= eip{(X)) = 
From ^3] we find 



2 
2 



Lp{x)dx, 0^ = 67/1(00) = - I ip{x)dx. (3.27) 



lim 



lim = 0, 



and at i = 00, 



(3.28) 



'Pi,a(oo) = 2c^, Qi,a(oo) = 2c^, 7?.i,a(oo) = 0, 
^i,a(oo) =Qia(cx)) = 7ti,A(cx)) = 0. 
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Hence 



(ai ) = We, (ai ,/32) = - c^, (ai,/33) = g£ 

(«2,/3l) 
( "2 , ^2 ) 
( "2 , ) 
( "3 , /?! ) 
( "3 , ^2 ) 
( "3 , ^3 ) 



[^1,A - ai,A (1 - Ve)] , 

2 A 

[7^l,A + ai,A (ge - C^)] , 

Z A 
z A 



■Pl,A - fiia (1 - We) 
7^1, A + ai,A (<?£ - C^) 
7^1, A + ai,A 



(3.29) 
(3.30) 

(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 



As an illustration, let us do the computation that led to H3.31|l in detail. As in 
|14| . we use the facts that S'* = S, and {S + Si?*) X ~ R which can be easily 
seen by writing R = X]fcli('5'^)''- Furthermore we write R{x,ak) to mean 

\im^Rix,y). 
ye. J 

In general, since all evaluations are done by taking the limits from within J, we 
can use the identity XSk — Sk inside the inner products. Thus 

(a2,/?2) - Y [{{S + RS) {1 - XX) ,St - Soo) + a,.x ( , <5t - <5oo )] 
A 

[{{1-XX),{S + R'S) {5t -6oo))+ ai.A {Qe{t) - Q.(oo))] 

A 

= J [{{l-XX),{S + R'S)X {5t-5^)) +ai,x{qe-c^)] 
A 

= ^ [( (1 - A X) , Rix, t) - R{x, oo) ) + ai,A (q^ - c^)] 
A 

= Y [T^i.^it) - 7?.i,a(oo) + ai^A (qe ~ c,p)] 
A 

= Y [T^i.x{t) + ai^x {qe - c^)\ ■ 

A 



Wc want the limit of the determinant 

dot (^j/c-(aj,/?/c))i<j-fe<i, 



(3.36) 



as — > oo. In order to get our hands on the limits of the individual terms 
involved in the determinant, we will find differential equations for them first as 

in [m. 
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Row operation on the matrix show that ai^x and di^\ fall out of the deter- 
minant; to see this add A 01^^/(2 A) times row 1 to row 2 and Xdi^x/X times 
row 1 to row 3. So we will not need to find differential equations for them. Our 
determinant is 



det 



XVi, 



X 



\ 



2 A 



1 



2 A 



2 A 



(3.37) 



A A 'A 

Proceeding as in jjl] we find the following differential equations 



d 

^Qr,A = .,.(A-l-7^,.), 



d 

= g„ - g„ «e -Pn 

^"^l.A =P« (A-7^l,A), 
^^1,A = 2l,A - 9« A, A, 

^'Pl,A = p„ (^A - 1 - 7^1, A 



(3.38) 
(3.39) 
(3.40) 
(3.41) 



Let us derive the first equation in H3.39|l for example. From ^21 (equation 2.17), 
we have 

— = -R{x,t) 



Therefore 

dQi,x 
dt 



d 

It 



Q{x,t)dx ~ [1- X) I Q{x,t)dx 



— 00 
t 



dQ 
dt 



dx-{l- A) 



dQ 

at 



dx 



i?(a;,t) da; - (1 - A) g„ + (1 - A) g„ / R{x,t)dx 



(1 - X)R{x,t)dx 



Now we change variable from t to s where t = r(s) = 2aVN+ j^^- Then we 
take the limit N 00, denoting the limits of qe,Vi,\, Qi,a, 7?.i,a, "Pi, a, Qi,a, "^i.a 
and the common limit of We and Ve respectively by g, Pi^a, Qi,a, 7?.i,a, ^i,a, Qi,a, 7?.i,a 
and u. We ehminate Qi^x and Qi^x by using the facts that Qi,a = ^i,a + A V2 
and Qi.A — 'Pi, a- These limits hold uniformly for bounded s so wc can inter- 
change lim and -j^. Also limTv^oo N~^^^Qn — limw^oo N~^/^p^ = q , where q 
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is as in p.9|l . We obtain the systems 



(3.43) 

^Pi,A = g (1 - A - ^1,;,) , = -9 (3.44) 

The change of variables q fi — q{x) dx transforms these systems into 

^^r,. ^ ^ (^r,. - A) , A^,. ^ -L (2^,. + , (3.46) 

^Pi.A = --^ (l - A - 1i,a) , = V2fi,x- (3.47) 

Since hmg^oo = 0, corresponding to the boundary values at i = oo which 
we found earlier for Pi.a, ^i,a, ^i.A, "^i.A, we now have initial values at fi — 0. 

Therefore 

Pi,a(0) = ^i,a(0) = ^i.a(O) = ^i,a(0) = 0. (3.48) 

We use this to solve the systems and get 

^=^e^ + ^e- (3.49) 
2^/2 2^/2 

..^^e--^.-.^. (3,50) 

— - A , ^1+ A 

ni,x = ^ — ^ 2 — ^ (3-^2) 

^i,A = ^ (^'^ - ^'')^ ^i^A = ^ + - 2)- (3.53) 

Substituting these expressions into the determinant gives H2.1|l . namely 

A - 1 -cosh/i(s,A) + \/a sinh/x(s,A) 
Z;i(s, A) = Z;2(s, A) ^—^ , (3.54) 

where Df} ~ limAr^oo Di3,n- 
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3.3 The GSE Case 

The GSE case is the easy one. All calculations in and go through 
essentially unchanged except for the trailing factor of A. Therefore we will not 
reproduce them here. 

3.4 Interlacing property 

The following series of lemmas establish Corollary (|2.2|) : 
Lemma 3.3. Define 



I A 

a., = 



^ dX3 \ 2-X 
Then Uj satisfies the following recursion 



(3.55) 



A=l 



1 */ j = o, 

U - 1) for j > I, j even, (3.56) 

j Qj-i for j >l, j odd. 



Proof. Consider the expansion of the generating function /(A) = y around 
A= 1 

Since aj = jlbj, the statement of the lemma reduces to proving the following 
recurrence for the bj 

f 1 if J = 0, 

bj = < ^ for i > 1, j even, (3.57) 

[_ bj-i for j > 1, j odd. 

Let 



These are the even and odd parts of / relative to the reflection A — 1 —(A — 1) 
or A ^ 2 — A. Recurrence (|3.57f) is equivalent to 

A.f^-^"(x)^^x-i)A.fo<id(^x) 

which is easily shown to be true. □ 
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Lemma 3.4. Define 



/(s,A) = l 



2- A 



tanh 



/i(s, A) 



(3.58) 



for X = 2X- A^. Then 



^2 n 



9A2 



A=l 



I g2n+l 

2n+ 1 aA2"+i 



1 i/ 71 = 0, 

ifn>l. 



(3.59) 



Proof. The case n = is readily checked. The main ingredient for the general 
case is Faa di Bruno's formula 



dt 



(hit)) 



]_dh\ ' 



1 d''h 
n\ dP^ 



, (3.60) 



where k = X]r=i ^^'^ above sum is over all partitions of ?i, that is all 
values of fci , . . . , A:„ such that X]r=i "^^i = n. We apply Faa di Bruno's formula to 
derivatives of the function tanh tLh^^ which we treat as some function g{X{X)). 

Notice that for j > 1, \x-i nonzero only when j = 2, in which case it 
equals —2. Hence, in (|3.60() . the only term that survives is the one corresponding 
to the partition all of whose parts equal 2. Thus we have 

92"-fc ^(s,X) 
tanh ^^ ' 



dX 



2n-k 



A=l 



dx^ tanii^— 



if fc = 2j + 1, j > 
for k = 2j, J > 



^(.s,A) 
tanh 



Qy2n+l-k 



A=l 







(_l)i-J (2n+l-k)\ 



tanh^ 



if = 2j, j > 
for k = 2j + 1, j > 



A=l 



Therefore, recalling the definition of aj in 13.55|l and setting k = 2j, we obtain 

2n 



:2n 



9A2 



■/(s,A) 



A=l 



2ri\ a* 



fc / (9A'= V 2- A dX^-^-'^ 



tanh ^^ ' ' 



E 

fc=0 

2^ -TTT^TW-^ a2.. tanh ■ 



A=l 



j=0 



(2j)!(n-j)! dX^-3 



A=l 
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Similarly, using fc = 2 j + 1 instead yields 



i2n+l 



/(s,A) 



A=l 



2n+l 

E 

k=0 



2n+l\ d 



tanh ■ 



k J dX'' \ 2-X aA2"+i-'= 
(2j)!(^-j)! 2j + l aA"-J- 



tanh 



2 

^(s, A) 



A=l 



since a^-j^^/{2i + 1) = a^.. Rearranging this last equality leads to H3.59() . □ 
Lemma 3.5. Let Di(s,A) and Di{s,\) be as in H2.1|l and H2.2() . Then 

fi{s, A) \ 



Di{s,X) = D4s,\) 1- 



A 

2- A 



tanh ■ 



(3.61) 



Proof. Using the facts that —1 — cosh a; = — 2 cosh^ f ' ^ ~ cosh^ x — sinh^ x and 
sinh X = sinh | cosh f we get 



„ / , ^ —2 ~, ~,A+vA sinh /i(s, A) 

Di{s,X) = ——D4{s,X) + D2is,X)— '^^ ' ' 



A(s,A)+Z?2(s,A) 



A-2 

-2 
A-2 

^4(S,A) 

-D4(s,A) 
Di{s,X) ( 1-2 



A-2 



A cosh^ + A sinh^ + ^A sinh /x(s, A) 



A-2 



1)4(5, A) A sinh^ + Va sinh/x(s. A) 



cosh 



2 /' mI^.a) 



A-2 



A sinh^ ''^ 2"^' + 2 sinh 




cosh(^(^/)) 


(A - 2) cosh^ 




) 



2- A 



tanh' 



fi{s,X) X 



^4(5, A) tanh^ 



2 2-A 
2 



tanh' 



^i(s. A) 



□ 



For notational convenience, define di(s. A) = dJ^^(s, A), ^4(5, A) = 1)^^(5, A). 
Then 



Lemma 3.6. For n > 0, 

Y Q2n+1 Q2n 



(2 71+1)! aA2»+i (2n)! ^A^' 



(ii(s, A) 



(-1)" 9" 



A=l 



i! 9 A' 



■ ^4(5, A) 



A=l 
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Proof. Let 



A , n(s, A) 

tanh ^ „ ^ 



2- A 



by the previous lemma, we need to show that 
2 g2n+i 2 92 ?i 



(2 71+1)! aA2»+i (2n)! ^A^ 
Now formula (|3.60() applied to ^4(5, A) gives 



d4(s,A)/(s,A) 



; ^ «4(s, A) 

n! 9 A" 



Qk 

Therefore 
1 



^4(5, A) 



A=l 







if fc = 2.7 + 1, j > 0, 



d4(s,A) iffc = 2j,i>0. 



(^;^a^'^4(s,A)/(,s,A) 



2n+l 



(2n+l)! 



fe=0 

(-ip 



A=l 



^= — dd „ , „ — ^ - , ^ f 

^ 0\9n — 9 1-1-1 



j^^{2n-2j + iy.j\ dX^ *aA2"-2j+i 



A=l 



Similarly 



1 ^2 



(2n)! 5A2' 



Therefore 
1 



1 ^" /2n\ d'' Q2n-k 

.d4(s,A)/(.,A) 

(2n-2j)!j! Jxi 



A=l 



j=0 



A=l 



2n+l 



1 



(2n+ 1)! 9A2"+i (2n)! SA^ 
(-ip 



E 



(2n-2j)!j! 



(i4 (s, A) 



d4(s,A)/(s,A) 

g2n-2j 



A=l 



^2n-2j+l 



QX27i-2j-' 2ri-2j + l 9A2"-2j+ 



Now Lemma 13.41 shows that the square bracket inside the summation is zero 
unless j = 71, in which case it is 1. The result follows. □ 



A=l 



Lemma establishes the inductive step in the proof of CoroUarv 12.21 
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4 Numerics 



Let 



qn{x) 



gn 



q{x, A) 



(4.1) 



SO that qo equals q from H1.7|) . In order to compute Fi3{s,m) it is crucial to 
know qn accurately. Asymptotic expansions for (7„ at — oo are given in ^3]. We 
outline how to compute qo and qi as an illustration. From |lHj . we know that, 
as t ^ +00 



90 

gi(-i/2) 



10657 13912277 



2t9 



8t 



12 



cxp 



:it3/2) 



17 



1513 



850193 



1 

407117521 



24t3/2 2^3H^ 2i"34t9/2 2i535i6 



O 



/ 1 



(4.2) 



Quantities needed to compute Fp{s,m),m 
integrals involving q^, such as 



1,2, are not only q^ but also 



{x — s) q^^{x) d X ^ Jo ~ qo{x)dx 



(4.3) 



Instead of computing these integrals afterwards, it is better to include them 
as variables in a system together with q^, as suggested in j^. Therefore all 
quantities needed arc computed in one step, greatly reducing errors, and taking 
full advantage of the powerful numerical tools in MATLAB. Since 



^0 = - 



ql{x)dx, III = ql, J'q = -go, 



(4.4) 



the system closes, and can be concisely written 



d_ 

ds 



( 10 \ 

lo 
V -^0 / 



sqo + 2^0 



V 



11^ 
1l 

-90 



(4.5) 



We first use the MATLAB built-in Runge-Kutta based ODE solver ode45 to 
obtain a first approximation to the solution of (|4.5|l between x = 6, and x = —8, 
with an initial values obtained using the Airy function on the right hand side. 
Note that it is not possible to extend the range to the left due to the high 
instability of the solution a little after —8; (This is where the transition region 
between the three different regimes in the so-called "connection problem" lies. 
We circumvent this limitation by patching up our solution with the asymptotic 
expansion to the left of x = —8.). The approximation obtained is then used as a 
trial solution in the MATLAB boundary value problem solver bvp4c, resulting 
in an accurate solution vector between a; = 6 and x = —10. 



19 





100 X 100 


100 X 400 


Fi-Percentile 


Ai 


A2 


A3 


Ai 


A2 


A3 


0.01 


0.008 


0.005 


0.004 


0.008 


0.006 


0.004 


0.05 


0.042 


0.033 


0.025 


0.042 


0.037 


0.032 


0.10 


0.090 


0.073 


0.059 


0.088 


0.081 


0.066 


0.30 


0.294 


0.268 


0.235 


0.283 


0.267 


0.254 


0.50 


0.497 


0.477 


0.440 


0.485 


0.471 


0.455 


0.70 


0.699 


0.690 


0.659 


0.685 


0.679 


0.669 


0.90 


0.902 


0.891 


0.901 


0.898 


0.894 


0.884 


0.95 


0.951 


0.948 


0.950 


0.947 


0.950 


0.941 


0.99 


0.992 


0.991 


0.991 


0.989 


0.991 


0.989 



Table 1: Percentile comparison of Fi vs. empirical distributions for 100 x 100 
and 100 x 400 Wishart matrices with identity covariance. 



Similarly, if we define 



{x ~ s)qo{x)qi{x)dx, Ji 



qo(x)qi(x)dx, (4.6) 



then we have the first-order system 



/ 91 \ 



d 
ds 



h 



J 



V 



9i 

s qi + &ql qi 

qoqi 

-90 91 



(4.7) 



which can be implemented using bvp4c together with a "seed" solution obtained 
in the same way as for go ■ Work is in progress to provide publicly downloadable 
versions of the MATLAB routines. 

Table ^ shows a comparison of percentiles of the Fi distribution to corre- 
sponding percentiles of empirical Wishart distributions. Here Ai denotes the i^^ 
largest eigenvalue in the Wishart Ensemble. The percentiles in the A^ columns 
were obtained by finding the ordinates corresponding to the Fi -percentiles listed 
in the first column, and computing the proportion of eigenvalues lying to the 
left of that ordinate in the empirical distributions for the A^. The bold en- 
tries correspond to the levels of confidence most commonly used in statistical 
applications. The reader should compare this table to a similar one in [Jj. 
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